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. Abstract: We analyze the state space of a Bianchi I universe with anisotropic sources. Here 

we consider an extended state space which includes null geodesies in this background. The 
evolution equations for all the state observables are derived. Dynamical system approach is 
used to study the evolution of these equations. The asymptotic stable fixed points for all the 
evolution equations are found. We also check our analytic results with numerical analysis of 
these dynamical equations. The evolution of the state observables are studied both in cosmic 
^ time and using a dimensionless time variable. Finally the cosmic microwave background 

anisotropy maps are generated, assuming that the universe is anisotropic and dominated by 
I one of the anisotropic sources since decoupling. We find that they contribute dominantly to 

, t>J[> CMB quadrupole. 

^ . 1 Introduction 

l/^ . In standard cosmology it is assumed that the space-time is homogeneous and isotropic. Af- 

ter the discovery of temperature anisotropics of the Cosmic Microwave Background (CMB) 
and accelerating expansion of the universe, a standard cosmological model describing universe 
dominated by cold dark matter (CDM) and cosmological constant (A) is formulated, known 
I as ACDm model or cosmic concordance model. Fluctuations in the temperature of CMB ra- 

diation are statistically isotropic in this model [HE]- A central assumption in our standard 
model of cosmology(and in most cosmological models) is that the universe is homogeneous and 
^ . isotropic upto small perturbations, and thus described by a perturbed Freidmann-Robertson- 

\^ ' Walker(FRW) metric. Obersevations show that the temperature of the CMB is isotropic to a 

^ ■ remarkable degree indicating that our universe is close to an Friedmann-Lemaitre (FL) model. 

According to Ellis, Geren and Sachs (EGS) theorem [3], if the CMB temperature were exactly 
isotropic about every point in spacetime, the universe would have to be exactly an FL model. 
This result is not directly applicable to CMB, because, CMB radiation is not exactly isotropic. 
The presence of this feature in the observations of temperature anisotropics of the CMB those 
from Wilkinson Microwave Anisotropy Probe(WMAP) seems to be inconsistent with a homo- 
geneous and isotropic FRW model. It may be evidence that we live in a globally anisotropic 
universe. 
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The current CMB data supports an inflationary Big Bang model of cosmic origin for our 
universe. However at large angular scales, there are some anomalies observed in CMB data, such 
as, a low value for quadrupole power, alignment of quadrupole and octopole modes roughly in the 
direction of virgo cluster, ecliptic north-south power asymmetry, parity asymmetry between even 
and odd multipoles, almost zero correlations on large angular scales of CMB and an anomalous 
cold spot ( 10° diameter) in the southern galactic hemisphere of CMB sky [3]. Several other 
solutions to low multipole CMB anomalies can be found, for example in Ref. |5j. Some these 
supposed deviations seen in CMB data were addressed by the WMAP science team in Ref. [6]. 

Generally, in a global anisotropic universe, during inflation shear decreases and eventually 
it goes over to an isotropic phase with negligible shear. In order to produce any substantial 
amount of shear in recent times one needs to induce anisotropy in the spacetime. One way to 
induce anisotropy is to have anisotropic matter present at the last scattering surface. It was 
shown earlier that power suppression in CMB quandrupole without affecting higher multipoles, 
can be explained by assuming anisotropic matter [71[5]. The sources for these matters can be 
uniform magnetic field, cosmic strings or a domain walls. A nanogauss scale magnetic field could 
be present today which would have been produced during inflation due to a lorentz-violating 
term to the photon sector [9]. We call matter in this form as "Lorentz- Violating magnetic 
field (LVMF)". 

In this paper our approach is not to deal with metric approach directly, but rather to adhere 
to an approach similar to orthonormal frame formalism [TO]. This paper is organised as follows. 
In section [21 we write down the Einstein's equations in cosmic time. Geodesic equations in 
cosmic time are studied in section [31 In section [H we do a fixed point analysis of all the 
evolution equations in terms of a dimensionless time variable (r), to study the asymptotic 
evolution of the state observables. Finally, in section [HI we show the temperature patterns for 
CMB due to various anisotropic sources considered in our study. 

2 Bianchi I universe 

The way we study dynamical systems is that we first write evolution equations in terms of 
few dimensionless state observables. It is a very invaluable tool to obtain qualitative informa- 
tion about the solution about the state space of Bianchi universes ((see ref. [TO] and references 
therein). Here we concentrate only on Bianchi 1 model. First we derive the Einstein equa- 
tions in metric approach. Then we write in terms of state observable analogous to orthonomal 
frame formalism( [11]. In this approach one writes the field equations as first order differential 
equations. An advantage of this approach is that the derivation of geodesic equation is easy 
compared to the metric approach. 

We start with a Bianchi I line element with a residual planar symmetry in the y — z plane 

as 

ds'^ = dt^ - aitfdx^ - h{tf{dy'^ + dz^). (1) 

We choose the diagonal energy momentum tensor of the form Tjf = (p, —pa, —Ph-, ~Pb)y where 
Pa = WaP and pb = WbP- The four kinds of anisotropic matter we study here are given in Table 

m 

Einstein equations for the metric given in Eq.[T] are 

2HaHb + Hl = SttGp, (2) 
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Matter 


Wa 


Wb 


Magnetic Field 


-1 


1 


Cosmic String 


-1 





Domain Walls 





-1 


LVMF 


1 






Table 1: Matter content for different anisotropic sources considered in this paper. 
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Figure 1: H' vs. t'. In the left figure, h'{0) = 173, ho/Ho = 0.01, JIq = 0.1,p'(0) = 10, and in 
the right figure h'{0) = -173,/io/i^o = 0.01, J^o = 0.1,/>'(0) = 10. 



2Hb + 3Hl = -STTGpa, 
Ha + Hb + Hl + HaHb + Hl = -STTGpb, 

with Ha = ^ ^ and Hb = j;^- The equation of continuity is given by, 

p+{Ha + 2Hb)p + HaPa + 2HbPb = 0. 



(3) 
(4) 



(5) 



Here now onwards we choose to set SttG = 1 for the rest of this paper. This can be otherwise 
be seen as rescaling, all the quanties in Eq. [2, 3, 4] with l/^/SnG. Now, we rewrite the above 
set of equations Eq.[2, 3, 4], [5] in terms of the average expansion H = ^°-~^^^b ^^^^ shear 



Ha— Ml 



'' . The Einstein's equations and the equation of continuity in cosmic time becomes 



^ = -H^--^h^--^{p + Pa + 2pb), 

dh , 1 , 

Tt = ^hH + -j={pb-pa), 

dp oui , Pa + 2pb. , o ^ / ^ 

Tt = -mp+^^)+2^{p,-Pa) 



(6) 
(7) 
(8) 



and the constraint equation can be expressed as 



(9) 
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Figure 2: h' vs. t' . In the left figure, /i'(0) = 173,/io/i?o = 0.01, = 0.1,p'(0) = 10, and in 
the right figure /i'(0) = -173,/io/i^o = 0.01, f^o = 0.1,/>'(0) = 10. 
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Figure 3: p' vs. t' . In left{right) figure, the initial data is same as in left{right) figures of Fig. 
[1], respectively. 
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We now introduce dimensionless variables H' , h' , t' defined as 



H' = ^, (10) 

t' = tHo, (12) 

P' = ^ (13) 
Po 

where Hq and /iq are the current value of hubble and shear paprameter respectively. Then the 
equations for the state observables in terms of these dimensionless variables can be written as 

H' = -H'^-l^ho/Hofh'^-^{l + Wa + 2wh)nop', (14) 

h' = 3h'H' + ^{wb-Wa)nop', (15) 

P = -3H{1 + )p +2—-j={wb-Wa)p , (16) 

where overdot represents differentiation with respect to t' and = irwi- Now, we solve the 
above equation numerically for all the cases listed in Table. [T] The evolution of the dynamic 
variables H',h' and p' are plotted in Fig. [1], [2] and [3]. For all these anisotropic matter 
sources, it turns out that H' decreases with increassing time t' . Uniform magnetic field and 
LVMF matter have almost similar rate of expansion, where as walls and strings have a higher 
rate of expansion. Similarly from Fig. [2] we see that isotropy can be attained at late times for 
all cases of anisotropic matter and shear becomes negligible at late times. Here we find that 
isotropy sets in much later for walls compared to the other three sources. Except for walls, in 
all other cases isotropy is achieved almost at the sametime. We observe from Fig. [3] that, the 
energy density also decreases much faster for magnetic field, LVMF and strings compared to 
walls. From Fig. [5], it can be observed that evolution of shear depends on the initial conditions. 
Signature change of shear can be seen for certain kind of anisotropic matter. This signature 
change may have implications in early universe. From the left plot of Fig. [2], we see that this 
change occurs for walls, whereas right plot of Fig. [2] shows a similar signature change but 
for cosmic strings, depending on a positive/negative shear in the beginning. However late time 
behaviours of H',h',p' in cosmic time are independent of initial conditions. 

3 Null Geodesies in cosmic time t' 

As we know the geodesic equations are of second order differential equations in metric approach. 
If we reperesent the same equation in terms of H and h, then the geodesies can be written as 
first order differential equations which are given by [12] 

^ - He - (17) 
e - -He - , (17) 

-H + 2^y,, (18) 
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k2 = [-H-^k2, (19) 
4 = {-H-\\k^, (20) 



where e = /cq is the energy of photon, /ci, k2 and A;3 are the three components of photon 
momentum. The variables e and ki (z=l,2,3) satisfies the constraint equation 

= kik\ (21) 

In terms of cosmic time t' the above equations become 

In order to solve the null geodesies Eq. [25], one needs to take care that the solutions should 
satisfy the constraint equation Eq. [21]. We solve the geodesic equations numerically for two 
particular cases. 

Case 1 : In this case only ki is taken to be nonzero and k2 and are zero. Initial values 
of A;i(0) = 2QeV and A:2(0) = ^^3(0) = are taken. The constraint equation Eq. j21j gives 
ki = ibA^o = ±e. Here we choose the positive value for kQ. The plot for kQ are shown in left of 
Fig. [3] for only initial values of A;o(0) = 20. Here we find that if initial value of k^ is positive 
then for entire evolution it remains positive. It's also true for negative initial vaues of feg. Other 
intial values are same as Fig. [1], [2] and [3]. We observe that k^ decreases slowly for magnetic 
field and faster for walls. This can be observed from Fig. [1] where H' decreases slowly for walls 
and faster for a magnetic field. Higher expansion rate means lower energy. Hence the above 
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Figure 5: Shown here are left ki vs. t' and right vs. t' . The initial data is same as in Fig 
[T] and A;i(0) = A;3(0) = 25.3. These are plotted for Case 2. 

analysis is consistent. 



Case 2: In this case A;i(0) = A;3(0) = 25.3. Other intial values are same as Case 1. The 
plots for ki and fcs are shown in Fig. [5]. In this case we also observe that feo decreases slowly 
for magnetic field and faster for walls, where as ki decreases faster in case of walls and slowly 
in case of walls. 



4 Fixed point analysis of the evolution equations 

In the previous section we studied the evolution of state observables in terms of a time variable 
anolagous to cosmic time. In order to study the asymptotic evolution of our Bianchi I universe 
with various anisotropic sources considered here, we resort to dynamic systems approach. So, 
we cast all the evolution equations in terms of dimesionless variables which give us a set of 
dynamical equations. The advantage of working in this approach is that we can show the 
cosmological dyanamics insensitive to the intial conditions. Here we define the dimensionless 
variables t, a as dr/dt = H and a = h/^/3H. These can be called as expansion normalised time 
and shear variables. Thus our Einstein's equations and the continuity eqaution are given by 



dH 

Ik 
da 

Ik 



-H{l + q), (26) 



-3(1 -a') 



^(1 - w)a - {wb - w) 



(27) 



and 



+ cr^ = 1 (28) 

where q = 2a^ + {1 + 3w)/2n, w = {wa + 2wb)/3 and n = p/3H'^. 

Next, we write the (null) geodesic equations in terms of dimensionless variables Ki = 
ki/ko{i = 1, 2, 3). They are given by 

^ = (l + s)e, (29) 
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Matter 


a 


(Kl,i^2,i^3) 


Magnetic Field 


1 


(±1,0,0) 


Cosmic String 


1/2 


(±1,0,0) 


Domain Walls 


-2/5 


(0,±1/V2,±1/V2) 


LVMF 


-1 


(0,±l/x/2,±l/^) 



Table 2: Stable fixed points of the whole state space. 



Ki = {s + 2a)Ki, (30) 
K2 = {s-a)K2, (31) 
= {s-a)K^, (32) 

where s = (1 — 3Ki)a. The variables Ki (i=l,2,3) satisfy the constraint equation 

KiK' = 1 . (33) 

Here we will analyze the fixed points of the full set of equations combining Einstein equations, 
equation of continuity and the geodesic equations [12]. Fixed points for a can be easily obtained 
from Eq. [27] which are ±1,-1 and 2{wb — w)/{l — w). All the stable fixed points are given 
in Table [2]. Same conclusion can be drawn from solving these equations numerically. The 
evolution of a and Ki and are plotted in Fig. [6] and [7]. For positive r, we can see that 
a evolves towards its stable fixed points while Ki and evolves towards their corresponding 
stable fixed points. This whole analysis is independent of the initial conditions. Exact solutions 
correspopnding to fixed points of a can be derived with the method as described in appendix 
of ref. [13]. 



5 Temperature Patterns 

Assuming that the universe is a Bianchi type I and dominated by an anisotropic matter, we 
find the constribution to the CMB anisotropics due to these sources. It has been shown earlier 
that the total quadrupole anisotropy in the presence of a uniform magnetic field, can be made 
small compared to that obtained with isotropic sources in standard ACDM model [7]. The 
tempearture of the cosmic microwave background as function of angular coordinates (9, cj)) on 
the celestial sphere is given by |14j . 

pTO—TLSS gSAo 

T{e, (^) = Tlss I , (34) 

y^e^^o ±(l-e6^o)cos2e 

where T^ss is the temperature of the cosmic microwave background at the surface of last 
scattering, tq and tlss a.re the values of r at today and at the time of last scattering surface 
and 

/•TO 

Ao= adr . (35) 

J TLSS 

Here we have taken tq = and t^ss = —7. It turns out that the quadrupole contribution 
is dominant for all the four cases of anistopic matter considered here. The CMB quadrupole 
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Figure 7: Left Ki vs. r , and Right vs. r. Here initial value a = 



.01 is taken. 



Figure 8: Temperature patterns -p^gj for Magnetic field {top left), Cosmic Strings {top right), 
Domain walls {bottom left), and LVMF {bottom right). Here we have taken (t{0) = 10~^. 

temperature patterns are given in Fig. [5]. We find that all other higher multipoles receive 
negligible contribution due to these sources. That is why here we only show temperature 
patterns for / = 2 only. We also observe from Fig. [S] that temperature maps for magnetic field 
and cosmic strings are of similar patterns and whereas temperature maps for domain walls and 
LVMF show similar patterns. This can be understood from the evolution of a given in Fig. [6]. 
The value of a stays positive for the entire evolution in case of magnetic field and cosmic string 
where as it becomes negative for the entire evolution in case of walls and LVMF. 

6 Conclusions 

Here we study the evolution of Bianchi-I universe in the presence of different types of anisotropic 
matters. First we see the evolution of the state space variables H' ,h' , p' ,ko,ki,k^ in cosmic 
time. Depending on the initial conditions we can have sign change in observable like shear. 
This may have some applications in early cosmology. Then we examine the evolution in terms 
of dimensionless variables like a, Ki , in dimensionless time r. In the r— frame we find the 
fixed points of the full set of evolution equations including geodesic equations. We also check 
the stability of the fixed points numerically. In r time the temperature of the CMB can be 
written in a very simple manner. Finally we generate the temperature patterns of the CMB 
for four different types of anisotropic matter considered here. We find that the contribution to 
CMB temperature patterns mainly coming from quadrupole. We also pointed out the difference 
between the CMB tempearture maps due to all the cases of anoisotropic matter. 

Here we have not constrained parameters like shear in these models from CMB data. It will 
be interesting to generalise our result to other Bianchi class A and B models. In case of other 
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Bianchi models the state space would be larger and we need to find different subspace of state 
space and their stable fixed points. It will also be interesting to study the CMB polarisation 
anisotropy in these models. 
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